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Abstract. For matrix analogues of embedded surfaces we define discrete cur- 
vatures and Euler characteristics, and a non-commutative Gauss— Bonnet the- 
orem is shown to follow. We derive simple expressions for the discrete Gauss 
curvature in terms of matrices representing the embedding coordinates, and 
provide a large class of explicit examples illustrating the new notions. 



1. Introduction 

A particular way of discretizing surfaces by replacing functions by matrices has for 
a long time been used in physics to obtain a quantum theory of surfaces (mem- 
branes) moving in Minkowski space, sweeping out 3-manifolds of vanishing mean 
curvature |Hop82| . The discretization, sometimes called "Matrix Regularization" , 
is of independent mathematical interest and contains many interesting structures. 
One of the main features of the correspondence between functions on the surfaces 
and matrices is that the Poisson bracket of two functions becomes the commutator 
of two matrices. This allows for an easy construction of discrete analogues of any 
expression involving products and Poisson brackets of functions. In a recent paper 
AHHIQ], the geometry of surfaces embedded in Riemannian manifolds has been 
expressed in terms of Poisson brackets of the embedding coordinates. Using these 
formulas, one can define discretizations of the Gaussian curvature and the Euler 
characteristic, and it is immediate to prove a discrete Gauss-Bonnet theorem (see 
Theorem EG)]). 

Apart from being interesting in their own right, these discrete concepts might 
also help to solve questions related to the regularization in the above mentioned 
Membrane Theory. For instance, solving the equations of motion in Membrane The- 
ory yields matrices corresponding to a discrete surface. As solutions corresponding 
to surfaces of arbitrary topology exist, one would like to be able to determine the 
geometry from the matrices in some way. In Theorem 13.111 we provide formulas 
for computing the discrete curvature and the discrete Euler characteristic given 
the matrix analogues of the embedding coordinates (which are the solutions to the 
equations of motion in Membrane Theory). Thus, in the limit of large matrices one 
may determine the Euler characteristic, and hence the topology, of the surface. 

2. Surface geometry and Poisson brackets 

Let us recall some of the results obtained in fAHHlO . Consider a surface E em- 
bedded in a Riemannian manifold M, of dimension m = 2 +p, via the coordinates 
x 1 (u 1 , it 2 ), . . . , x m (u , u 2 ), where it 1 ,?/ 2 are local coordinates on E. Furthermore, 
let n\(v}, it 2 ), . . . , n™(u 1 ,M 2 ) for A = 1, . . . ,p denote the components of p orthonor- 
mal vectors Na normal to the surface at each point. Indices k, I will run from 
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1 to 771 and indices a, b will run from 1 to 2. The metric of M is denoted by g^, 
the Christoffel symbols by f* fc and the covariant derivative by V. Regarded as a 
subspace of TM, the tangent space TE is spanned by the vectors e a — (d a x l )di. 

Letting /^(u 1 , u 2 ) be an arbitrary non- vanishing density on £, one defines a Pois- 
son bracket on C°°(£) by setting 

(2.1) {f,h} = -e ah (d a f){d b h), 

where e ab is antisymmetric with s 12 = 1 and d a = With this bracket we define 
the tensors 

(2.2) V" {.<■■.. r'\ 

(2.3) SJ = - p e ab (d a x l )(V b N A ) 3 = {x\n\} + {x\x k }t{ in l A , 

and one can also consider them as maps TM — > TM by lowering the second index 
with the ambient metric <?, i.e. 

(2.4) V(X) = V'n^X'il 

(2.5) SA(X)=S i J°g hj X*a i . 
With these definitions, one finds that 

(2-6) TrSi = (S A y(S A ) j = -^det(h A ,ab) 

P 

(2.7) Tr V 2 =V\P\ = -2-^, 

where g = det (g(e a , eb)) is the determinant of the induced metric on E, and h AMb 
is the second fundamental form corresponding to the normal vector N A . 
In the main part of this article we will make use of the following result: 

Theorem 2.1 ([AHH10 ). Let K denote the Gaussian curvature o/E. Then 

1 o 2 P 

(2.8) K = -g(R(e 1 ,e 2 )e 2 ,e 1 ) - |- £ TrS 2 A , 

9 9 A=i 

where R is the curvature tensor of M. 



3. Matrix regularizations and discrete curvatures 

In the following, we shall assume that E is a compact closed orientable surface. Let 
us first define what is meant by a "matrix regularization" , and then show some of 
its properties. 

Definition 3.1. Let Ni,N2,... be a strictly increasing sequence of positive inte- 
gers, let {T a } for a = 1, 2, . . . be linear maps from C°°(E, K) to hermitian N a x N a 
matrices and let H(N) be a real-valued strictly positive decreasing function such 
that limTv-i-oo Nh(N) < oo. Furthermore, let ui be a symplectic form on E and let 
{•, •} denote the Poisson bracket induced by w. If {T a } has the following properties 



DISCRETE CURVATURE AND THE GAUSS-BONNET THEOREM 



3 



for all f,h E C°°(E) 

(3.1) lim ||T<*(/)|| <oo, 

a— >oo 

(3.2) lim \\T a (fh) -T a {f)T a (h)\\ = 0, 



(3-3) 
(3-4) 



lim 



±.[T a (f),T a (h)]~T a ({f, h} 



= 0, 



lim 2nh a TiT a (f) = [ fu, 



where || • || denotes the operator norm and h a = h(N a ), then we call the pair (T a , K) 
a matrix regularization of (Y*,lu). 

Given local coordinates u , u on S, we write u) = p^u 1 ,u 2 )du 1 A du 2 , and it is easy 
to see that the induced Poisson bracket becomes 

1 



{f,h} 



(d a f)(dbh). 



Definition 3.2. If fx, fa, ■ . ■ is a sequence of matrices such that f a has dimension 
N a and if it holds that 



(3.5) 



lim 

Oi — VGO 



fa ~ T a {f) 



0, 



then we say that the sequence converges to the function f . 
Definition 3.3. A matrix regularization (T a ,h) is called unital if 



(3.6) 



lim 111 



r(i) 



Q. 



Remark 3.4. Although unital matrix regularizations seem natural, and all our ex- 
amples fall into this category, it is easy to construct examples of non-unital matrix 
regularizations. Namely, let (T a ,h) be a matrix regularization and consider the 
map T a defined by 

/ 



f a (f) = 



T a (f) 



0/ 



Then (T a , K) is a matrix regularization which is not unital, since 



lim 



f a (l)-t Na+1 



> 1. 



Proposition 3.5. Let (T a ,h) be a unital matrix regularization. Then 
(3.7) lim 2irN a h a = u). 

Proof. Let us use formula (|3.4p with / = 1. 

w= lim 27rft Q TrT Q (l) = lim 27rft Q Tr |"r Q (l) + t N - Ijv, 
= lim (2nh a N a + 2^ Tr(T Q (l) - 1^)) = lim 2i:h a N a 
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since 



lim |27rft Q Tr(T Q (l) - t Na )\ < lim 2j:h a N a ||T Q (1) - t Na \\ = 0, 

due to the fact that the matrix regularization is unital. □ 

Proposition 3.6. Let (T a ,h) be a matrix regularization of (E,lu) and let {f£} be 
a sequence of matrices converging to fk € C°°(E) for k — 1, ...,n. T/ien {ai/" + 
■ • • + a n /n} converges to a\f\ + ■ ■ ■ + a„/„ /or any oi, . . . , a n G K and 



lim 

a— >oo 



lim 

a— >oo 



/r ■••/« I < TT Um n ra ^) 



fe=l 



(3.8) 

(3.9) 
(3.10) 



Proof. The first statement about {ai/f + • • • a„/"} follows directly from the lin 
earity of T a . Let us prove (|3.8p by induction on n. Thus, we assume that 
holds and compute 



lim 2irh a Tr (/« ■••/")= / A " " " /. 



lim 

a— »oo 



Pet J*Q 

/l * * * Jn+l 



< lim 



fl ■ • • /» /«+! " T a (/ n+1 ) + T a (/„+i) 



< lim 



lim 



?a Pa ra rriot i r \\\ , ro Pq \\rnot f £ \ 

Jn+l ~ J Un+lJ \\ + Jl ' ' ' J„ 



n+l 



/f •••/„" l|r Q (/„ + i)||< n lim ll TQ (/fc)ll' 



fe=i 



To prove (|3.9p we again proceed by induction and assume that (|3.9p holds for any 
given n, and then compute 



lim 

a— f oo 



1 ' ' ' Jn+l _ 1 1/1 *" Jn+l] 



< lim ( /?•••/« 



< lim 



< 



/n+i - r Q (/„ + i)| | + 1 |/f • • • fZT a (f n+1 ) -r(/ r .. / n+1 )| I ) 

1 ••■/„ Jn+l Un+lJ + Jl ■■■Jn 1 Un+lJ 

T a (/l • • ■ /n+l) - T Q (/l • • • fn)T a (f n+1 ) + T a (/i ■ ■ ■ /„)T Q (/„ +1 )| |) 

lim ( /?•••/„" f%+i-T a (f n+1 ) +||T»(/ n+ i)|| .//•••/;: /"!/:•••. A J 



+ l|r Q (/l • • • fn)T a (f n+1 ) - T a (/i • • • / n+1 )|| ) 

Finally, we prove the trace formula: 
lim 2wh a Tr (/f • • • /«) 



= 0. 



lim 2nh, 

a—>oo 



TrT Q (/i • • • /„) + Tr (/f •••/«- T Q (/i • • • /„)) 
/l • • ■ /-.w, 
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since 



lim 

Oi— > oo 



h a Tr(f?---t-T a (h---fn)) 



< lim H a N a 



by formula (133)) . 



□ 



The above result allows one to easily construct sequences of matrices converging to 
any sum of products of functions and Poisson brackets. Namley, simply substitute 
for every factor in every term of the sum, a sequence converging to that function, 
where Poisson brackets of functions may be replaced by commutators of matrices. 
Proposition 13.61 then guarantees that the matrix sequence obtained in this way 
converges to the sum of the products of the corresponding functions. 

Proposition 3.7. Let (T a ,h) be a matrix regularization and assume that {/«} 
converges to f. Then {f\} converges to f. 

Proof. Due to the fact that ||^4|| = \ \A^\\ one sees that 



lim 

a— >oo 



fi-T a (f) 

since {f a } converges to / 



lim 

a— >oc 



lim 

a— >oo 



fa ' T a (f) 



0. 



□ 



If the matrix regularization is unital, one can relate the matrix sequence converging 
to the function 1//, to the inverse of a sequence converging to /. 

Proposition 3.8. Let (T a ,H) be a unital matrix regularization and assume that f 
is a nowhere vanishing function and that {f a } converges to f . If f^ 1 exists and 
\\fa~ 1 \\ * s uniformly bounded for all a, then {f^ 1 } converges to l/f. 



Proof. One calculates 



lim 



f- 1 - T a (l/f) 



< lim 



fa 1 ^N a - f a T a {l/f) 



= lim 

Oi—^OO 

< lim 
= 0, 



t 1 l Na -f a T a (l/f) + T a (l)-T a (l) 



fa 



T a (l) 



f a T a (l/f)-T a (l) 



since the matrix regularization is unital and H/q, 1 !! is assumed to be uniformly 
bounded. □ 

Recall that E is embedded in a m = 2 + p dimensional manifold M via the embed- 
ding coordinates x (u , it 2 ), . . . , x m {u l , u 2 ), and that p orthonormal normal vectors 
are given with components n\. By {X l a } and {N l A a } we will denote arbitrary 
sequences converging to x l and n\ respectively. Moreover, given the metric g^ 
and the Christoffel symbols f * fc of M, we let {Gij. a } and {Tj k a } denote sequences 
converging to g^ and T^ k respectively. To avoid excess of notation, we suppress 
the index a whenever all matrices are considered at a fixed (but arbitrary) a. 
In analogy with (|2.3p we define 

(3.11) (S A ){ = ±-{X\N%]G k , k + l[^,X ; ]ff;>XG fe ^, 
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and 

(3.12) tr s% = {(s A y k )\s A ^. 

Let g a b be the induced metric on E, and g its determinant. We set 

(3.13) 7 =^, 

P 

and denote by {7 a } an arbitrary sequence of invertible matrices converging to 7. 
By denning 



in 



it follows from (|2.7[) that 



(3-14) -^(nM = ^G) k [X\Xi][X\X l ]G, 



2 y fc; 1 2h 2 



converges to 7 2 . 

If the embedding space is K m , the above formulas reduce to 

ra 

(3-15) trS% = -h£ [X\iVi][^,iVi], 

1 1 m 

(3.16) ira^N^El 1 '^ 2 ' 

and in R 3 one obtains 

(3.17) frS 2 = ^ £ i" £ *''^ t )- 1 [Jf i .[^^]][y,[^^'l]f 1 1 



since 



(3.18) ^ = —4^,^}, 



1 

27" 

defines a unit normal vector to the surface (cp. |AHH10j , where (13.171) is also given 
for arbitrary codimension). 

We are now ready to define and present formulas for the discrete curvature in a 
matrix regularization of E. 

Definition 3.9. Let (T a ,ti) be a matrix regularization of (E, u>) and let K be 
the Gaussian curvature of E. A Discrete Curvature of E is a matrix sequence 
{Ki, K2, K3, . . .} converging to K, and a Discrete Euler Characteristic of E is a 
sequence {xi, X2, X3, ■ ■ •} such that lim \a = X- 

a— >oo 

From the classical Gauss-Bonnet theorem, it is immediate to derive a discrete ana- 
logue for matrix regularizations. 

Theorem 3.10. Let (T a , K) be a matrix regularization o/(E,w), and let {K\, K2, ■ ■ .} 
be a discrete curvature of E. Then the sequence Xi)X2j • • • defined by 



(3.19) XQ = fi Q Tr 

is a discrete Euler characteristic o/E. 



la K a 



DISCRETE CURVATURE AND THE GAUSS-BONNET THEOREM 



7 



Proof. To prove the statement, we compute lim a _>oo \ a and show that it is equal 
to x(S). Thus 



lim x a = lim —2-Kh a Tr 

ct — ^ oo a— foo Z7T 



and by using Proposition 13.61 we can write 

lim Xa = 77- / K^-lo = 77- / K^-pdudv = ^- [ Ky/gdudv = x(S), 
«->-oo 27T J s p 2ir p 2ir J E 

where the last equality is the classical Gauss-Bonnet theorem. □ 

Theorem 3.11. Let (T a ,h) be a unital matrix regularization of and let 

{Ki 2 } be a matrix sequence converging to g(R(ei, e-i)ei, e-n / g (the sectional curva- 
ture of TT, in M). Then the sequence K of matrices defined by 

(3.20) K = K 12 - \ j^m-^trSDr 1 

is a discrete curvature o/S. Thus, a discrete Euler characteristic is given by 

(3.21) x = hTv (ffKn) ~ f E Tr [r^r&A ■ 

A=l 

Proof. By using the way of constructing matrix sequences given through Proposi- 
tion [3751 the result follows immediately from Theorem 12. II □ 

Note that if p = Jg, then 7 = 1 which implies that one can choose j a = l_/v Q when 
the matrix regularization is unital. 

4. TWO SIMPLE EXAMPLES 

4.1. The round fuzzy sphere. For the sphere embedded in R 3 as 
(4.1) x = (2; 1 , x 2 , x 3 ) = (cos<^sin#, sin<£>sin0, cos#) 

with the induced metric 



(4-2) ( gab ) 



1 
sin 2 6 



it is well known that one can construct a matrix regularization from representations 
of su(2). Namely, let Si, S2, S3 be hermitian N x N matrices such that [S^ , S k ] = 
ie^tS 1 , (S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 = (N 2 - l)/4, and define 

(4.3) X 1 = ■ 2 S\ 

Then there exists a map TW (which can be defined through expansion in spherical 
harmonics) such that T^Xx 1 ) = X 1 and (T^, H = 2/y/N 2 - 1) is a unital matrix 
regularization of (S 2 ,^/gdd A d(p) |Hop82 . A unit normal of the sphere in M 3 is 
given by N £ TM. 3 with N = x l di, which gives N l = X 1 , and one can compute the 
discrete curvature as 

1 rn 

(4.4) ^ N = ~fi? ^ Tr[X\X>] 2 = t N 

i<j=l 
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which gives the discrete Euler characteristic 

2N 

(4.5) xn = HTrK N = hN ■■ 



y/N 2 - 1 ' 
converging to 2 as N — > oo. 

4.2. The fuzzy Clifford torus. The Clifford torus in S 3 can be regarded as 
embedded in R 4 through 

x = (x 1 , x 2 , x 3 , x 4 ) = ^p(cos (fx, smtpi, cos ^2, sin ^2), 
v2 

with the induced metric 

1 A 



(5ab) - 2 l j > 

and two orthonormal vectors, normal to the tangent plane of the surface in TR 4 , 
can be written as 

N± = x 1 d 1 + x 2 d 2 ± x 3 <9 3 ± x 4 d 4 . 

To construct a matrix regularization for the Clifford torus, one considers the N x N 
matrices g and h with non-zero elements 

g kk =uj k - 1 for k = 1,...,N 
hk,k+i = 1 for A; = 1, .. . , N- 1 
= 1, 

where to = exp(i2#) and 8 — ir/N. These matrices satisfy the relation hg = ujgh. 
The map T^ N ' is then defined on the Fourier modes 

y_ — imi^i+im 2 v32 

as 

rW(Y^) = u}i m ^ m2 g mi h m2 , 

and the pair (T^ N \ h = sin#) is a unital matrix regularization of the Clifford torus 
with respect to y/gdipi A d<f2 [FFZ891 |Hop 89 . Thus, using this map one finds that 

*> , + ft > 
*>'-*> 

= ±X 3 and iV| = ±X 4 . By a 



X 1 - 


Tix 1 ) 


= ^^(COS^!) = 


x 2 = 


T(x 2 ) 


= ^ T ( S ' m( Pl) = 


x 3 = 


T(x 3 ) 


= ^ T ( COS ^2) = 


x 4 = 




= ^ T (sin^ 2 ) = 


which implies that N± 


= x\ 


Nl = X 2 , Nl -- 



straightforward computation one obtains 

4 



\ [X\Xi] 2 =21 



h 2 
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and therefore 

4 4 

— £ [X\Nl\[X\Nl\ =-^J2 [X\Xi? 
and since [X\X 2 ] = [X 3 , X 4 ] = it follows that 

4 4 

^ £ [X\Nl][X\Nl] = [X\X] 2 = 

This implies that the discrete curvature vanishes, i.e. 

4 -, 4 



^ = i E [^.^i]^',^;] + i ^ [X\Nl][X\Nl_] =1-1 = 0, 
which immediately gives \n = 0. 



5. AXIALLY SYMMETRIC SURFACES IN R 3 

Recall the classical description of general axially symmetric surfaces: 
(5.1) af = (/(it) cosv, /(it) sin v,h(u)) 

'h'(u) cosu, h'(u) sinv, — f'(u)), 



which implies 

, , //' 2 + ft' 2 0\ , v ±:i (h'f"-h"f' u 

(Sab) - n f 2 (.«abj = 



/V ^<"V ^72 + / ,2 V -fh' 

where h a b are the components of the second fundamental form. The Euler charac- 
teristic can be computed as 



1 f f u + h'(h'f" -h"f) _ f 



(5.2) x =— Ky/g = - J-Ldu=- 

which is equal to zero for tori (due to periodicity) and equal to +2 for spherical 
surfaces (due to f'(u±) — Too)- 

While a general procedure for constructing matrix analogues of surfaces embed- 
ded in R 3 was obtained in [ABH+09bllABH+09al (cp. also IArn08bj ). let us restrict 
now to h(u) = u = z, hence describe the axially symmetric surface X as a level set, 
C = 0, of 

(5.3) C(x) = ±(x 2 +y 2 -f 2 (z)), 

to carry out the construction in detail, and make the resulting formulas explicit. 
Defining 

(5.4) {F(x),G(x)} R s = VC- (VFx VG), 
one has 

(5.5) {x,y} = -ff'(z), {y,z} = x, {z,x} = y, 
respectively 

(5.6) [X, Y] = ihff'{Z), [Y, Z] = ihX, [Z, X] = ihY 
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for the "quantized" ( "non-commutative" ) surface. In terms of the parametrization 
given in (|5.ip . the above Poisson bracket is equivalent to 

(5.7) {F(u, v), G(u, v)} = e ab (d a F) (d b G) 

where d\ — d v and di = d u . By finding matrices of increasing dimension satisfying 
(|5.6j) , one can construct a map T a having the properties (|3.2j) and (|3.3j) of a matrix 
regularization restricted to polynomial functions in x, y, z (cp. [Arn08a ). 

For the round 2-sphere, f(z) — 1 — z 2 , ()5.6j) gives the Lie algebra su(2), and its 
celebrated irreducible representations satisfy 

(5.8) X 2 + Y 2 + Z 2 = 1 if h 



\/N 2 - 1 



When / is arbitrary, one can still find finite dimensional representations of (|5.6j) as 
follows: rewrite (15.61) as 



(5.9) [Z, W] = hW 

(5.10) [W,W^] = -2hff(Z) 

implying that z^ — Zj — h whenever Wij ^ and Z diagonal. Assuming W = X+iY 
with non-zero matrix elements Wk^k+i = Wk for k = 1, . . . , N — 1, one thus obtains 
(with wq = wyv = 0) 

Z kk = ^(N+l-2k) 

wl ' = -2hff(h(N + 1 - 2k)/2) = Q k , 

which implies that 

i=i 

and the only non-trivial problem is to find the analogue of (|5.8p . To this end, define 

(5.11) f 2 = X 2 + Y 2 = -(W f + W*W), 

with W given as above. As Z has pairwise different eigenvalues, the diagonal matrix 
given in (|5.11[) can be thought of as a function of Z; hence as f 2 {Z). It then trivially 
holds that 

(5.12) C = X 2 + Y 2 -f 2 (Z) = 0, 

for the representation defined above. The quantization of h comes through the 
requirement that f 2 should correspond to f 2 . While for the round 2-sphere f 2 
equals f 2 , provided h is chosen as in ()5.8|) . it is easy to see that in general they can 
not coincide, as 

[X 2 + Y 2 ~f(Z) 2 , W] = \{WW ] + W ] W)/2 - f(Z) 2 , W] 

= 1 W [ W \W] + ^[W\W]W - f(Z)[f(Z), W] - lf(Z), W]f(Z) 

= ■■■ = f(z) {hf(z)w - [f(z),w]) + (hf(z)w - [f(z), w])f(z) 

with off-diagonal elements 
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that are in general non-zero (hence X 2 + Y 2 + f 2 (Z) is usually not even a Casimir, 
except in leading order). 

How it does work is perhaps best illustrated by a non-trivial example, f(z) = 
1-z 4 : 

h 4 / 

(5.13) w 2 k = — [(N + lfk - 3(N + l) 2 k(k + 1) + 

2(N + l)k(k + l)(2k + 1) - 2k 2 {k + l) 2 ) 

fk = \U + = ^((N + If (2k 1) - 6(N + l) 2 k 2 

+ A(N + l)k(2k 2 + 1) - 4fc 2 (fc 2 + 1)) 
(note that Wq = w 2 N = is explicit in (|5.13l0 so that 

4 (AT J- 1^3 



(5.14) iX' + Y' + Z' ).. =K 



Expressing the last two terms via Z 2 (note that the cancellation of fc 3 and k 4 terms 
shows the absence of Z 3 and higher corrections) one finds 

X 2 + Y 2 + Z 4 + h 2 Z 2 = h 4 (iV + 1)2 f (N + l) 2 - 4(N + 1) + 4)1 

16 V / 

16 



which equals 1 if h is chosen as 2/ \/N 2 — 1. Note that this is the same expression 
for H then for the round sphere, f 2 = 1 — z 2 (cp. (|5.8p ). 

A more elegant way to derive the quantum Casimir (cp. also |Roc91| [GPS09 ) 

(5.15) Q = X 2 + Y 2 + Z 4 + h 2 Z 2 

is to calculate 

[X 2 + Y 2 + Z 4 , W] = [(WW* + W^W)/2 + Z 4 , W] 
= --- = h 2 [W,Z 2 }, 

which determines the terms proportional to H in the Casimir. 
Due to the general formula 

(5-16) K = -±- iEjueim tf)-*[X\ [X\X 1 ]} [X*, [XV,X«]}r 2 

one obtains, for the axially symmetric surfaces discussed above, 

(5.17) K = r 2 {[.f.f'?(Z) + ±_[W, ff(Z)]Wl + ±Wi[W, ff\Z)]^r 2 
with 

(5.18) f = + W*W) + (ff') 2 (Z) = f(Z) 2 (f'(Z) 2 + 1) + 0(H), 
giving 

(5.19) k = -(f'(Z) 2 + ty 2 f(z)- 1 f"(z) + 0(h) 
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and for f(z) 2 = 1 — z 4 one has 

(5.20) K = (AZ 6 + t- Z 4 ) ~ 2 (6Z 2 - 2Z 6 ) + 0(h) 

(5.21) 7 2 = 1 - Z 4 + AZ 6 + 0{H). 

Note that (cp. (|5.9|0 zj — Zj-i — h for arbitrary /, and that (due to the axial 
symmetry) K and j 2 are diagonal matrices, so that 



in this case simply being a Riemann sum approximation of J K^/g, indeed converges 
to 2, the Euler characteristic of spherical surfaces. 
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